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^ . I. INTRODUCTION 

P.. 

The determination of the properties of hadronic matter at finite temperature and den- 
sity is a fundamental problem in nuclear physics. In particular, the study of liquid-gas 
phase transition in medium energy heavy-ion collisions is of considerable interest. Many 
intermediate-energy collision experiments have been performed [1] to investigate the un- 
known features of the highly excited or hot nuclei formed in collisions [2,3]. Theoretically, 
much effort has been devoted to studying the equation of state for nuclear matter and to 
discussing the critical temperature, T^. The calculated critical temperature of symmetric 
nuclear matter lies in the range 13 - 24 MeV for various phenomenological models [4]- [12]. 
Glendenning [13] first discussed the phase transition with more than one conserved charge 
and applied the method to the possible transition to quark matter in the core of a neu- 
tron star. Miiller and Serot [14] discussed asymmetric nuclear matter using the stability 
conditions on the free energy, conservation laws and the Gibbs criterion for the liquid-gas 
phase transition. The liquid-gas phase transition of asymmetric nuclear matter has also 
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been discussed in effective chiral models [15,16]. It was found that the critical temperature 
decreases with increasing asymmetry parameter, a. 

For finite nuclei, there is another temperature which is called the limiting temperature, 
as pointed out by Levit and Bonche [5]. Below the limiting temperature, nuclei can 
exist in equilibrium with the surrounding vapor. When the temperature is higher than TJj^ 
the nuclei arc unstable and will fragment. This is called Coulomb instability. The size cff^ect 
and Coulomb interaction are important in the determination of the limiting temperature, 
resulting in a lower limiting temperature compared with that for infinite nuclear matter [17]. 
Recently, Natowitz et al. obtained the limiting temperature by using a number of different 
experimental measurements [18]. Prom these observations the authors extracted the critical 
temperature of infinite nuclear matter Tc = 16.6±0.86 MeV [19]. Their results show that the 
limiting temperature is in good agreement with the previous calculations employing either 
a chiral symmetric model [20] or the Gogny interaction [8]. 

To study the properties of hadronic matter, we need phenomenological models since 
QCD cannot yet be used directly. The symmetries of QCD can be used to determine largely 
how the hadrons interact with each other. On this basis, models based on SU{2)l x SU{2)r 
symmetry and scale invariance were proposed. These effective models have been widely used 
in recent years to investigate nuclear matter and finite nuclei, both at zero temperature and 
at finite temperature [20]- [22]. Papazoglou et al. extended the chiral effective models to 
SU{3)l X SU{3)r including the baryon octets [23,24]. As well as the models based on the 
hadron degrees of freedom, there are additional models based on quark degrees of freedom, 
such as the quark meson coupling model [25,26], the cloudy bag model [27], the NJL model 
[28] and the quark mean field model [29]. Recently, we proposed a chiral SU{3) quark mean 
field model and applied to investigate hadronic matter and quark matter [30]- [33]. This 
model is very successful in describing the properties of nuclear matter [30], strange matter 
[31,32], finite nuclei and hypcrnuclei [33] at zero temperature. A successful model should 
describe well the properties of nuclear matter, not only at zero temperature but also at finite 
temperature. In this paper, we will apply the chiral SU (3) quark mean field model to study 
the liquid-gas phase transition and Comloub instability of asymmetric nuclear system and 
compare our results with the recent experimental analysis. 

The paper is organized as follows. The model is introduced in section II. In section III, 
we apply it to investigate nuclear matter at finite temperature. The numerical results are 
discussed in section IV and section V summarises the main results. 

II. THE MODEL 

Our considerations are based on the chiral SU{3) quark mean field model (for details 
see Refs. [31,33]), which contains quarks and mesons as basic degrees of freedom. Quarks 
are confined into baryons by an effective potential. The quark meson interaction and meson 
self-interaction are based on SU{3) chiral symmetry. Through the mechanism of sponta- 
neous chiral symmetry breaking, the resulting constituent quarks and mesons (except for 
the pseudoscalars) obtain masses. The introduction of an explicit symmetry breaking term 
in the meson self-interaction generates the masses of the pseudoscalar mesons which satisfy 
the partially conserved axial- vector current (PCAC) relation. The explicit symmetry break- 
ing term of the quark meson interaction leads in turn to reasonable hyperon potentials in 
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hadronic matter. For completeness, we introduce the main concepts of the model in this 
section. 

In the chiral limit, the quark field q can be split into left and right-handed parts and 
1r- Q = Ql + Qr- Under SU (3) l x SU (3)/j they transform as 



LqL, q'u = RqR. 
The spin-0 mesons are written in the compact form 



(1) 



M(M+) = E ± m = ^ ^ (ct" ± ct") A", 

V2 a=0 



(2) 



where a" and vr" are the nonets of scalar and pseudoscalar mesons, respectively. A" (a = 
1, 8) are the Gell-Mann matrices, and A° = ^ I. The alternatives, plus and minus signs 
correspond to M and M+. Under chiral SU{3) transformations, M and transform as 
M ^ M' = LMR^ and M+ M+' = RM+L'^. In a similar way, the spin-1 mesons are 
introduced through: 



(3) 



with the transformation properties: /, 



RvnR'^. The matrices S, 



n, Vn and yl^ can be written in a form where the physical states are explicit. For the scalar 
and vector nonets, we have the expressions 



v2 „=o 



K* 



*o 



K 



*o 



(4) 



Pseudoscalar and pseudovector nonet mesons can be written in a similar fashion. 
The total effective Lagrangian has the form: 

^eS — ^qO + ^qM + -^ES + ^VV + ^xSB + -^Am^ + '^fe,+'^c, 



(5) 



(6) 



where Cqo = q ij'^d^ q is the free part for massless quarks. The quark-meson interaction CgM 
can be written in a chiral SU (3) invariant way as 

^qM — Qs 

(8 8 \ / S 8 \ 

E ^aXa + i E ^aXal'] ^ ' E T^^^^^ - E 7^7'<AJ ^. (7) 

a=0 a=0 J 2V2 Va=0 a=0 / 

In the mean field approximation, the chiral-invariant scalar meson >Css and vector meson 
Cvv self-interaction terms are written as [31,33] 
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= -\ kox' {a' + e) + k, (a' + + k, + C^) + k,x<y\ 



4 Xo 3 f^oCo 



2xg 



(8) 
(9) 



where 5 = 6/33; (Tq, Co and xo ai'e the vacuum expectation values of the corresponding mean 
fields (7, ( and x- The Lagrangian C^sb generates nonvanishing masses for the pseudoscalar 
mesons 



-XSB 



r 

Xo 



(10) 



leading to a nonvanishing divergence of the axial currents which in turn satisfy the relevant 
PCAC relations for tt and K mesons. Pseudoscalar and scalar mesons as well as the dilaton 
field, Xi obtain mass terms by spontaneous breaking of chiral symmetry in the Lagrangian 
(8). The masses of the d and s quarks are generated by the vacuum expectation values of 
the two scalar mesons a and C,. To obtain the correct constituent mass of the strange quark, 
an additional mass term has to be added: 



(11) 



where S — o (7 - As VS) = diag(0, 0, 1) is the strangeness quark matrix. Based on these 
mechanisms, the quark constituent masses are finally given by 



rrid 



V2 



and 



-9sCo + Am^, 



(12) 



where gs and Arris are chosen to yield the constituent quark mass in vacuum - we use 
m„ = rrid = 313 McV and = 490 MeV. In order to obtain reasonable hyperon potentials 
in hadronic matter, it has been found necessary to include an additional coupling between 
strange quarks and the scalar mesons a and ( [31]: 



Ch = {ha + /i2 C) ss . 



(13) 



In the quark mean field model, quarks are confined in baryons by the Lagrangian Cc = 
—^Xc^ (with Xc given in Eq. (14), below). The Dirac equation for a quark field under 
the additional infiuence of the meson mean fields is given by 



-ia-V + Xcir) + /3m* 



(14) 



where a = 77, /3 = 7 , the subscripts i and j denote the quark i {i = u, d, s) in a baryon 
of type j {j = N, A, S, S) ; Xci^^) is a confinement potential, i.e. a static potential providing 
confinement of quarks by meson mean-field configurations. The quark effective mass, m*, 
and energy, e*, are defined as 



iO 



(15) 
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and 



e* = ei-5>-5;0-5>, (16) 

where is the energy of the quark under the influence of the meson mean flelds. Here 
rriiQ = for i = u,d (nonstrange quark) and = Arris = 29 MeV for i = s (strange 
quark). Using the solution of the Dirac equation (14) for the quark energy e* it has been 
common to define the effective mass of the baryon j through the ansatz: 

M; = ^Ef- < p*l^ > , (17) 

where E* = + Ejspin is the baryon energy and < pfcm > is the subtraction of the 

contribution to the total energy associated with spurious center of mass motion. In the 
expression for the baryon energy Uij is the number of quarks with flavor "i" in a baryon 
with flavor j, with j = N {p, n} , S {S^, S°} , S S^} , A and Ej^pm is the correction to 
the baryon energy which is determined from a flt to the data for baryon masses. 

There is an alternative way to remove the spurious c. m. motion and determine the 
effective baryon masses. In Ref. [34], the removal of the spurious c. m. motion for three 
quarks moving in a confining, relativistic oscillator potential was studied in some detail. It 
was found that when an external scalar potential was applied, the effective mass obtained 
from the interaction Lagrangian could be written as 

M;^J2r^,je*-El (18) 

i 

where E^ was calculated to be only very weakly dependent on the external field strength. 
We therefore use Eq. (18), with Ej a constant, independent of the density, which is adjusted 
to give a best flt to the free baryon masses. 

Using the square root ansatz for the effective baryon mass, Eq. (17), the conflning po- 
tential Xc is chosen as a combination of scalar (S) and scalar-vector (SV) potentials as in 
Ref. [31]: 

Xc{r)^^-[xl{r)+xt''{r)] (19) 

with 

xUr)^\Kr\ (20) 

and 

xf(r) = ifc,r2(l + 7°). (21) 

On the other hand, using the linear definition of effective baryon mass, Eq. (18), the confining 
potential Xc is chosen to be the purely scalar potential Xc(^)- The coupling kc is taken as 
kc—1 (GeV fm~^), which yields root-mean-square baryon charge radii (in the absence of a 
pion cloud [35]) around 0.6 fm. 
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III. NUCLEAR MATTER AT FINITE TEMPERATURE 



Based on the previously defined quark mean field model the Lagrangian density for 
nuclear matter is written as 

-S'c^V'TmV'^'' - gpi^BlnTzil^ff + Cm, (22) 

where 

Fi,u = d^uj,, - d^ujf, and p^^ = - (23) 

The term Cm represents the interaction between mesons which includes the scalar meson 
self- interaction >Css, the vector meson self-interaction Cyv and the explicit chiral symmetry 
breaking term C^sbi all defined previously. The Lagrangian includes the scalar mesons cr, C, 
and X) and the vector mesons uj and p. The interactions between quarks and scalar mesons 
result in the effective nucleon mass Mj^, The interactions between quarks and vector mesons 
generate the nucleon- vector meson interaction terms of equation (22). The corresponding 
vector coupling constants and Qp are baryon dependent and satisfy the SU (3) relationship: 

At finite temperature and density, the thermodynamic potential is defined as 



N=p,n 



where £'^(A;) = \J M'^^ + 't\ The quantity z/jv is related to the usual chemical potential, 
fiN, by z/jv = Pn ~ q!^^ ~ 9^ P- The energy per unit volume and the pressure of the system 
are respectively e = VL — + ujsiPn and p = — f2, where piv is the baryon density. 

The mean field equation for meson 0j is obtained by the formula dVL/d4)i — 0. For 
example, the equations for cr, C are deduced as: 

^ ^ 3(7 Xo 

( xV 2'9m^ ( xV 2dmp dM^ - 



kox'C - 4^1 {a' + C') C - ^hC' - k,x<y' - + ^ {^rnlF, - ^m^F. j = (26) 
where 

1 M* 

<^^>^— dk—-^ [nn{k) + nn{k) + np{k) + np{k)] . (27) 
TT^ Jo rj*[k) 

In the above equation, nq{k) and nq{k) are the nucleon and antinucleon distributions, re- 
spectively, expressed as 
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""^^^^ " exp [{E*{k) - vg) /ksT] + 1 ^^^^ 

and 

^g(^) = i(T?*(i\ / \ /I T^i I 1 = (29) 

exp [[E*[k) + Vq) /kBT\ + 1 

The equations for the vector mesons, uj and p, are: 

+ Ag^u^ + 12^4a;p2 = {pp + p„) , (30) 

Xo 

— mjp + 4c/4P^ + 12^40; V = o9!^{Pp - Pn), (31) 
Xo -J 

where Pp and p„ are the proton and neutron densities, expressed as 

Pq^^ dkk^ [nq{k) - nq{k)] {q=p,n). (32) 

TT^ Jo 

In order to describe asymmetric nuclear matter, one can introduce the asymmetry parameter 
a which is defined as 

a = (33) 

Pn 

where pn = Pn + Pp- For symmetric matter a = 0, while for neutron matter a = 1. 

Let us now discuss the liquid-gas phase transition. For asymmetric nuclear matter we 
follow the thermodynamic approach of Refs. [13] and [14]. The system will be stable against 
separation into two phases if the free energy of a single phase is lower than the free energy 
in all two-phase configurations. This requirement can be formulated as [14] 

F{T, p) < (1 - X)F{T, p') + AF(r, /), (34) 

with 

p = (1 - A)p' + Ap", < A < 1, (35) 

where F is the Helmholtz free energy per unit volume. The two phases are denoted by a 
prime and a double prime. The stability condition implies the following set of inequalities: 

If one of the stability conditions is violated, a system with two phases is energetically favored. 
The phase coexistence is governed by the Gibbs conditions: 

p;(T,p') = p;'(r,p"), {q-n,p), (38) 

p'(T,p')=p"(T,p"), (39) 
where the temperature is the same in the two phases. 
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IV. NUMERICAL RESULTS AND DISCUSSIONS 



A. Liquid-gas phase transition 



The parameters in this model are determined by the meson masses in vacuum and the 
properties of nuclear matter which were listed in table I of Ref. [36]. We first discuss the 
liquid-gas phase transition of symmetric nuclear matter. In Fig. 1, we show the pressure 
of the system versus nucleon density at different temperatures using the square root ansatz 
for the effective nucleon mass (Eq. (17)). At low temperature, the pressure first increases 
and then decreases with increasing density. The p — Pn isotherms exhibit the form of two 
phase coexistence, with an unphysical region for each. At temperature T = 15.82 MeV, 
there appears a point of inflection, where dp/dp^ = 0, d'^p/dp% = 0. This temperature is 
called the critical temperature. Symmetric nuclear matter can only be in gas phase above 
this temperature. The pressure versus nucleon density with the linear definition of effective 
nucleon mass (Eq. (18)) is shown in Fig. 2. In this case, the critical temperature is 17.9 MeV. 
In both cases, the calculated critical temperature is close to the recent result = 16.6 ±0.86 
MeV, which was obtained by Natowitz et al. [19]. The critical temperature calculated with 
the Walecka model is larger than 20 MeV. This large Tc is partially caused by the large 
incompressibility modulus K (~ 540 MeV). For the same reason the stiff Skyrme interaction 
SKI gives a large which is close to 20 MeV, while the soft Skyrme interaction SKM* gives 
a small Tc [7] . The QMC model and the effective model based on SU (2) chiral symmetry 
provide reasonable values of around 15-16 MeV [15,37]. 




For the asymmetric case, the situation is more complicated. One cannot get the crit- 
ical temperature from the p — Pn isotherms. The chemical potentials of the proton and 
neutron are different. We show the chemical potential versus asymmetric parameter a at 
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temperature T ~ 10 MeV and pressure p — 0.12 MeVfm"^ with the square root definition 
of effective nucleon mass in Fig. 3 (For convenience, we use the reduced chemical potential 
which is defined as JIn = Hn — Mjy.) The solid and dashed lines are for proton and neu- 
tron respectively. The Gibbs equations (38) and (39) for phase equilibrium demand equal 
pressure and chemical potentials for two phases with different concentrations. The desired 
solution can be found by means of the geometrical construction shown in Fig. 3, which 
guarantees the same pressure and chemical potentials for protons and neutrons in the two 
phases with different asymmetry parameter a. 
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FIG. 2. The pressure of symmetric nuclear matter p versus nucleon density pjv at different 
temperatures with the linear definition of effective nucleon mass. 

The pairs of solutions found using the method just described, yield a binodal curve which 
is shown in Fig. 4. There is a critical point where the pressure is about 0.205 MeV-fm~^ and 
the corresponding asymmetry parameter is around 0.7. The two phases have the same a 
and therefore the same density at this point. The binodal curve is divided into two branches 
by the critical point. One branch corresponds to the high density (liquid) phase, the other 
corresponds to the low density (gas) phase. Assume the system is initially prepared in the 
low density (gas) phase with a = 0.55. When the pressure increases to some value, the 
two-phase region is encountered at point A and a liquid phase at B with a low a begins to 
emerge. As the system is compressed, the gas phase evolves from point A to C, while the 
liquid phase evolves from B to D. If the pressure of the system continues to increase, the 
system will leave the two- phase region at point D. The gas phase disappears and the system 
is entirely in the liquid phase. This kind of phase transition is different from the normal first 
order phase transition where the pressure remains constant during phase transition. If the 
initial asymmetric nucleon gas is larger than some value, the system enters and leaves the 
two phase region on the same branch. For example, the system becomes unstable at point 
A' and a hquid phase with a higher nucleon density begins to emerge at B'. The system is 
compressed at a fixed total a, with the gas phase evolving from A' to C and the liquid phase 
from B' to D'.The system leaves the two-phase region point C which is still in the original 
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gas phase. Therefore, for a given temperature, if the total asymmetry parameter of the 
system is larger than a critical value, the system cannot change completely into the liquid 
phase, however large the pressure. In other words, for a system with a fixed asymmetric 
parameter a, there exists a critical temperature, above which the system can only be in the 
gas phase at any pressure. 



> 




FIG. 3. Geometrical construction used to obtain the chemical potentials and asymmetric pa- 
rameters in the two-phase coexistence at temperature T = 10 MeV and p = 0.12 MeVfm"^. 
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FIG. 4. Binodal curve at temperature T = 10 MeV. The points A through D and A' through D' 
denote two kinds of phase transition. 

The a dependence of the critical temperature is shown in Fig. 5. The solid and dashed 
lines correspond to the square root ansatz and the linear definition of effective nucleon mass. 
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respectively. Tc decreases with increasing a. For the square root case, when a is less than 
0.2, the decrease of is very small. When a is larger than 0.6, decreases very fast. If a is 
larger than 0.88, the system can only be in the gas phase at any temperature. In the linear 
case, the critical temperature is 2 MeV larger than that in the square root case. The hquid- 
gas phase transition can occur for nuclear matter with any a if the temperature is lower 
than the critical temperature. The critical temperature of asymmetric nuclear matter was 
also studied in the Walecka model, the derivative scalar coupling model and the QMC model 
[10,37] where the authors found different critical temperatures for protons and neutrons. The 
lower critical temperature was chosen to be the of the system as an approximation. Their 
results show that Tc decreases almost linearly with increasing a, which is different from the 
results shown in Fig. 5. We use the stability conditions on the free energy, conservation laws 
and the Gibbs criterion for the liquid-gas phase transition of asymmetric nuclear matter, 
following Miiller and Serot [14]. This is expected to be better than the approximation used 
in the earlier discussion. 
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FIG. 5. The critical temperature Tc versus asymmetric parameter a. The solid and dashed lines 
are for the square root and linear cases, respectively. 

B. Coulomb instability 

In this subsection, we discuss the liquid-gas phase equilibrium for finite nuclei. Compared 
with the case of infinite nuclear matter, the size effect and Coulomb interaction are important 
for finite nuclei. When the Coulomb interaction is considered, the chemical potential for the 
protons will have an additional term corresponding to a uniformly charged sphere: 

I^Coul =5-^' W 

where Z and R are the charge number and radius of finite nuclei. Meanwhile, the pressure 
also has an extra term: 
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Pcouiip) = ^^P, (41) 

where A is the mass number of the nucleus and p is its density. For a hquid droplet, the 
surface pressure is 

Psurf{T,p) = -2-f{T)/R, (42) 
In the above equation, j{T) is the surface tension suggested by Goodman et al. [38] 



7(T) = (1.14MeVfm-2) 



ST 



T 



(43) 



where Tc is the critical temperature for infinite symmetric nuclear matter. The Gibbs con- 
ditions for two phase equilibrium now becomes 

p(T, PL, ai) + PcouiipL) + PsurfiT, pl) = p{T, pv, av), (44) 

Pn{T, PL, oil) = Pn{T, pV , Oiy), (45) 

Pp{T, PL, ul) + Pcoui{pl) = Pp{T, pv, av). (46) 

When the temperature is higher than a temperature (limiting temperature), the above 
equations have no solution. The finite nuclei cannot exist in equilibrium with the surrounding 
vapor. 

We show in Fig. 6 the mass number dependence of the limiting temperature, T^j^, for 
nuclei along the line of /3-stability: 

Z = 0.5A - 0.3 X IQ-^A^/^. (47) 

The solid and dashed lines correspond to the square root ansatz and the linear definition of 
the effective nucleon mass, respectively. The experimental values obtained recently by Na- 
towitz [18] are also plotted in the figure for comparison. The calculated results in this model 
are in good agreement with the experimental data. The limiting temperature decreases with 
increasing mass number. This means that when the temperature is higher than the limiting 
temperature, the heavy nuclei will fragment to light nuclei. Two useful parameterizations of 
Tiim/Tc, valid for 10 < ^ < 208, are {Tu^/Tc) = 0.611 - 0.00193^ + 3.32 x 10"^^ (square 
root case) and {Tum/Tc) = 0.591 - 0.00203A + 3.80 x lO'M^ (linear case), which are com- 
parable with that given in Ref. [19]. Calculations from other models also show that Tum 
decreases with increasing mass number. Numerical results show that there is a relationship 
between T;j^ and T^.. For larger T^, the calculated is also larger. Therefore, the results 
for Tiim found with the Walecka model and the SKI interaction are larger, while the results 
of the SKM* interaction are smaller compared with the experiments [37]. The values of Tum 
calculated for both the QMC model and the effective model suggested by Furnstahl et al. 
[20,37] are close to our results. 
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FIG. 6. The limiting temperature T;j^ versus mass number A of finite nuclei. The solid and 
dashed lines are for the square root and linear cases, respectively. The points with error bars are 
from Ref. [19]. The data derived from the double isotope yield ratio and thermal bremsstrahlung 
measurements are represented by the filled circles and open squares, respectively. 

V. SUMMARY 

In this paper, we extended the chiral 5*^7(3) quark mean field model to finite temper- 
ature and density. This model describes properties of infinite nuclear matter, finite nuclei 
and hypernuclei very well at zero temperature. The saturation properties and compression 
modulus of nuclear matter are reasonable. The hyperon potentials arc close to the empirical 
values for hadronic matter. The results for finite nuclei and hypernuclei are also consistent 
with experiment. We therefore want to know whether this model can also describe the sys- 
tem at finite temperature. The liquid-gas phase transition of infinite nuclear matter and the 
Coulomb instability of finite nuclei at finite temperature are discussed in this model. All 
the parameters have been determined in earlier papers and there is no further parameter 
to be adjusted. The critical temperatures for two-phase coexistence of symmetric nuclear 
matter, Tc, is 15.82 MeV (square root case) and 17.9 MeV (linear case). Both of these values 
are close to the recent experimental value, 16.6 ± 0.86 MeV. Tc is found to decrease with 
either increasing asymmetry parameter, a, or increasing mass number for finite nuclei. The 
critical temperature, Tc, in the linear case is about 2 MeV larger than that in the square 
root case. The values in both cases are in good agreement with the experimental values 
found in Refs. [18,19]. 
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